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We consider the generahzed set of theories of gravitation whose La- 
grangians contain the term B? : L = y/—g{R + j3R^). Inserting the RW 
metric with an imposed non-singular and inflationary behaviour of the scale 
factor a(t), and using a arbitrary perfect fluid, we study the properties of p 
and p in this context. By requiring the positivity of the energy density , as 
well as real and finite velocity of sound, we can obtain the range of values of 
(3 that ensure the inflationary behaviour and absence of singularity. 

I. INTRODUCTION 

There are many motivations to study quadratic-order curvature Lagrangians in cosmol- 
ogy: quantum corrections and renormalization of divergences, [|]-|^ effective actions of 
superstring theory, inflation without scalar fields, [0 etc. 

The aim of the present work is to explore the consequences of adding quadratic-order 
curvature terms to the Einstein action. From the beginning, a specified non-singular and 
inflationary scale factor and the Robertson- Walker metric are imposed, and the unknown 
energy-momentum tensor is then analysed. 

The only energy-momentum tensor T^^ choice in this work will be an arbitrary perfect 
fluid described by the energy density p and the isotropic pressure p. Just a couple of simple 
energy conditions will be used here. 
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p>0, {p + p)>0, (1) 

and additionally the adiabatic speed of sound, = {dp/ dp), must obey 

0<vl<l. (2) 

For each value of the spatial curvature k, the form of p and p are obtained as functions of 
the scale factor a{t) and other constants. The energy (|l|) and velocity of sound @ conditions 
can give constraints on the constant f3 associated with the quadratic-order curvature terms. 
It is also possible to compare the violations of the conditions with and without the 
quadratic-order curvature terms, which gives some insight on the consequences of using 
quadratic Lagrangians. 

II. THE MODIFIED FRIEDMANN EQUATION 

Considering the fact |1[ that the quantity ^/^{R^.upaR^'''"^ - ^R^.yR^" + i?^), known 
as Gauss-Bonnet term, is a total divergence in four-dimensional spacetime, the R^^p^R^'^^'^ 
term can be expressed in terms of R^ and R^^R^^ , so the Riemann tensor is not necessary 
to account for the curvature square corrections. 

Representing a general quadratic order curvature Lagrangian (but without the cosmo- 
logical constant), let us take the Lagrangian 

L = v^(i? + PiR^ + (32Rp.R'"') . (3) 

Calculating the variation of this action with respect to the metric g'^", the modified Einstein 
equations are 

Gpu + PiW^^'J + = KTp,. (4) 

We will use the Robertson- Walker (RW) metric 



ds^ = -df + a{tf 



(5) 



where the spatial curvature k takes the values 0, +1,-1 (corresponding to flat, closed and 
open spatial sections). In the case of this spherically symmetric homogeneous metric, the 
Weyl tensor has null components. So the identity C^yp^C^^'P'^ = R^upaR^"'"' ~'2.RpuR^'" + \R^ 
is equal to zero, yielding another equation that now allows us to write Rf^^R^'^ in terms of 
: 



-g{R + l3R^). (6) 



With the use of a perfect fluid, the modifled Friedmann third-order differential equation 
is obtained : 

a'{tf + k Kp{t) 18/3 



a{tf 3 a{tf 



[k"" -?>a'{tf -a{tfa"{tf- (7) 



2 a'itf [k - ait) a" it)] + 2 a{tf a\t) a"'{t)] = , 

where /? = (Pi + y) different from zero means that the quadratic-order curvature corrections 
are present. 

With the RW metric (^), the energy-momentum conservation equations yield 

P W = > (8) 

ct 

so that p{a) can be calculated using p{a). 



III. NON-SINGULAR SCALE FACTORS WITH INFLATIONARY 

BEHAVIOURS 

Non-singular and inflationary scale factors for the three values of spatial curvature k 
are studied. They show inflationary behaviour matching the classical radiation dominated 
behaviour. 

I) The A; = case : 
The following equation M 
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allows to obtain a scale factor a(t), which for a = (or a » a^) yields a{t) = a^y/t 
(the known solution for radiation dominated era in FRW cosmology with A; = 0), and for 
t — oo (or a « ao) gives a{t) = aoexp(t/2a). Figure |T] shows this behaviour compared 
to a{t) = ao\/t. 

a(t) 




t 



FIG. 1. a{t) X t for the case k = 0, with uq = a = 1. The solid line is the non-singular scale 
factor with inflationary behaviour; the dashed line is the known scale factor, a{t) = ao\/t, for 
radiation era in FRW cosmology. 



It is possible to derive Eq. (^ in respect to t and write 

^ a(t) gp^ [ajty + ap^] 
^ ' 2[a(t)4 + a(t)2ao2 + ao4a] ' ^ ' 

deriving this Eq. (|10|), the second derivative can be obtained, 

II u\ _ ^(^) «o^[a(t)^ + ao^][3a(t)^ + qq^] « - a{tf a^'^ [a(t)^ + a^^f 

4[a(t)4 + a(t)2ao2 + ao^a]' ' ^ ' 

and deriving this Eq. (|TT]), the third derivative is given by 



a'" it) = ^(t) Ht)' + a,'] r 4_ 
^' 8[a{ty + a{tyao^ + ao^af^ v;[v; oj 

2 a{ty ao'^ [a{ty + Oq^] [ 15 a{t)^ + 15 a{tf a^^ + 4 Oq^] a + 

ao®[15a(t)^ + 12a(t)2ao^ + aoV^} • 
The scale factor a{t) and its derivatives are smooth functions of t. 

II) The k = +l case : 

Generahzing the equation Eq. (j^), 



(12) 



1-^1-^lto + alog 



a{tf 



a{tY + ao^ 



(13) 



allows to obtain a(t), which for a = (or a ^ ag) yields a(t) = aoyl — (1 — t/t^Y (the 
known solution for radiation dominated era in FRW cosmology with k = +1), and for 
t — oo (or a << ao) gives a(t) = aoexp(t/2a). 

As in the /c = 0, it is possible to derive Eq. ([T3|) in respect to t and write 



a'{t) 



a(t)ao \J — ciity [aity + 



(14) 



a{t)HQ + a{ty aoHo + 2aQ^^aQ^ - a(tf a 

and the second and third derivatives can also be obtained without problems. All these 
functions are smooth. 

Ill) The k = -I case : 



Also generalizing the equation Eq. 



-1 + Jl + 



to + a log 



a{tf 



(15) 



allows to obtain a(t), which for a = (or a » a^) yields a{t) 



-1 + (1 + t/tc 



(the known solution for radiation dominated era in FRW cosmology with k = —1), and for 
t — oo (or a << ao) gives ait) = ao exp(t/2a;). 

As in the /c = 0, it is possible to derive Eq. (p!5| ) in respect to t and write 
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a'{t) = a(t)ao [ajt)^ + ao^/^ ^^^^ 

and the second and third derivatives can also be easily obtained. All these derivatives are 
smooth functions. 



IV. ANALYSIS OF p AND P 

Instead of trying to obtain solutions for a{t) of the modified Friedmann equation, Eq. 
(0), which can only be solved numerically, another way of analysis will be chosen. Replacing 
the derivatives of a(t), then the energy density p{t) can be written as p(a), also depending 
on the constants k and (3 and the constants of the definition of a {t). 

I) The A; = case : 

The energy density is given by 

, , 3 ao^ (a^ + ao^)^ 27a^ oq^^ (a^ + ao^)^ a/3 

p{a) = 2 1 6~ ^ (-'-') 

4 (a^ + + cio^ a) k 2 (a^ + + oq^ C() >^ 

14 + 29 + 2 oq^ (10 - 3 a) - 5 (-1 + a) 

and a thorough analysis shows that p > if a and (3 are inside the region 

Qi = {0 < a < 1, /3 > U {a > 1, /?_ < /5 < , (18) 

where /?_ and /?+ are functions of a. The fact of adding the non-linear terms in the La- 
grangian is represented by /? 7^ 0, and does not force p to be negative. 

The Figure || shows the behaviour of p(a), which has no divergences if a > 0. 
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0.5 1 1.5 2 

FIG. 2. p{a) X a for the case k = 0, with uq = k = 1, a = 0.5, /? = 0.01. The sohd hne uses the 
non-singular scale factor with inflationary behaviour; the dashed line is obtained using the known 
scale factor, a{t) = ao\/t, for radiation era in FRW cosmology. 



The pressure p can be obtained from the Eq. (|), 

, , ao^ (a^ + ao^) [«® + 2 + Oq^ (1 - 7 «) - 3 
»(a) = o 

■f^*^' x//l,09, /1\0 



+ 



(19) 



9a^ao^^ (a^ + ao^)a /3 . ^12 + 459 ^10 2 g ^8 4 _ gg + 
2(a4 + a2ao2 + ao^a)^/€^ 

OCT _ 12 / 1 , „ \ „ _6 _ 6 / o^7 , 1 1 r, „ N , 



25 ao'^ (-1 + a) a - 6 a'' ao^ (-87 + 119 a) + 
3 ao® [ 70 + a (-207 + 22 a)] + ao^° [35 + a (-220 + 87 a)\ } , 

If a > 0, the pressure p is always negative for some region in the a(t) or t domain, 
independent of /3. For a > 0, p has no divergences. 

Testing whether p + p > is true just restricts a and (5 inside the region 



g2 = {0<a<l,/3_</5</?+} , 

where /3_ and /5+ are another functions of a, different from the analysis of p. 
The velocity of sound, 



(20) 



V 



dp dp /da 



(21) 



dp dp I da 

does not diverge if a and /3 are inside the same region Q2- But it can be imaginary for any 
value of (a, /3) 7^ (0,0). 

So the conclusion for the = case is that the addition of high-order terms in the 
Einstein action does not lead to any other anomaly in the effective perfect fluid, i.e., with 
or without /5, the behaviour oi p, p, p + p and Vs is the same (under the condition that a 
and [3 are inside the region Q2). 

II) The k = +1 case : 

Analogous to the A; = case, but with more difficult calculations, the energy density is 
given by 

p{a) = _ _ rQ6^2a"ao^ + a^ao"+ (22) 



[a^ + a^ + 2 a ag^ V — + Aq^ Q^) 
4 \l —cfi + cio^ (o-^ + cto^) a — 4 oq^ + 4 



+ 



216 ao^ a /? 



fi (a^ (a - ao) ao^ (a + ao) x 
(a^ + a^ a^^ + 2 oq^ yj—a^ + Oq^ a) k 

(a^ + ao^)^ (51 + 14 oq^ - 13 a^ qq^ + 8 ao^) a + 4 ao® (a^ - a ao^)^ x 
(-17 + 2ao^) - 16ao^° {-a^ + ao^)^a^ + \l -o? + ao^ x 
(a^ - ao^) (a^ + ao^)^ (s a*' - 34 - 33 a^ - 10 ao*') + 
8a^ao^ (7a^-2ao^) (a^ + ao^)^ + 48 ao® (-a^ + ao^) | , 

and a long analysis shows that p > if a and (3 are inside the region 

Qi = {0 < a < ttma.,/?- < /? < 0} U {a > a„,ax,/3 < 0} , (23) 

where (3- is function of a and ao, and Omax is a function of ao, ftmax ~ 0.73 aO. Figure |^ 
shows the region Qi. The non-linear terms in the Lagrangian do not force p to be negative. 
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FIG. 3. The region Qi for the case k = +1, where p > 0, with ao = 1. 



The pressure p can be obtained from the Eq. (P) and will not be shown here. Even p > 
is possible if a and f3 are inside the region 

Q2 = {a > an,in, (3 <f3-Up> p+} , (24) 

where /3_ and /3+ are functions of a and oq, and a^i^i is a function of Oq, amin ~ 0.84 aO. As 
/? = is not included in the region Q2, this result is a benefit of adding non-linear terms in 
the Lagrangian. 

To obtain p + p > 0, a and P must be inside the region 

Q3 = {0 < a < a^ax, /?-</?< 0} U {a > a„,ax, /5 < 0} , (25) 

where /?_ and /3+ are another functions of a, different from the analysis of p, and Omax is a 
function of oq, Omax ~ 0.80 aO. The regions Qi and are similar. 

The velocity of sound does not diverge if a and P are inside the same region Qs- But it 
can be imaginary for any value of [a, (3) 7^ (0,0); maybe only in the Q2 region Vg has real 
values, this must be confirmed. 

So the conclusion for the k = +1 case is that the addition of high-order terms in the 
Einstein action does not lead to any other anomaly in the effective perfect fluid, i.e., with 
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or without /3, the behaviour of p, p, p + p and Vg is the same (under the condition that a 
and P are inside the region Qs). 

There is even an advantage because p can be always positive, with only the problem that 
the region Q2 leads to almost suppression of the radiation era. The calculus of the region 
Q2 is not easy, and it is a task for further work to verify if inflationary behaviour can exist 
with p > 0, p > and < Vs < 1. 

Ill) The case : 

Like the k — +1 case, the calculations are very long, and the energy density is given by 

3 oo^ (a^ + ao^) [o^ - 4 Qq^ o? + c? f Oq^ - 4 ^Jo? + oq^ ol) 
p{a) = ^ ^ ^ + (26) 



216 ao^ a /3 f e 2 /.i 2 , 00 2\ ^ 

a qq (51 a + 32 oo ) x 



I 5 .0 c / n n\ 4 



(a^ + ao^) a - 4 oq^ (a^ + oq^) (17 + 2 ao^) - 16 (a^ + ao^) x 
ao^°«^ + (a^ + ao^)^^^ a^(a^ + ao^)^ (5 + 17a^ oq^ + 10 ao"^) - 
Sa^ao^ (a^ + ao^) (70^ + 200^) - 48 oq^ a^] } , 

and p > is impossible for a > 0, independent of {3. 

The pressure p will not be shown here. Like p, p > is impossible for a > and any (5. 

Obviously, for a > and any /3 it is not possible to obtain p + p > 0. 

The velocity of sound is divergent and imaginary if a > 0, independent of (3. 

The case k — —1 does not provide useful conclusions, because even without the high- 
order terms in the Einstein action there are a lot of anomalies in the effective perfect fluid. 
The hope that quadratic-order terms could avoid these anomalies was not achieved. 

V. CONCLUSION 

In all three cases of spatial curvature, it is possible to choose the values of a and (3 so 
that the quadratic-order terms in the Lagrangian do not add more anomalies to the effective 
perfect fluid. 
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For closed spatial section, the opposite happens, i.e., anomalies are suppressed, and this 
case of non-singular inflationary scale factor with well-behaved perfect fluid deserves more 
investigation. 

Other interesting improvements would be : use another solutions of the scale factor, 
more energy conditions, compare the (3's constraints with observational values and study 
many types of effective perfect fluids. 
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